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I. INTRODUCTION 

Exploitation of biological resources as practiced in fishery, and forestry has strong impact on dynamic 
evolution of biological population. The over exploitation of resources may lead to extinction of species which 
adversely affects the ecosystem. However, reasonable and controlled harvesting is beneficial from economical 
and ecological point of view. The research on harvesting in predator-prey systems has been of interest to 
economists, ecologists and natural resource management for some time now. 

The optimal management of renewable resources has been extensively studied by many authors [1,2, 
3, 7, 8, 12]. The mathematical aspects of management of renewable resources have been discussed by [10]. He 
had investigated the optimum harvesting of logistically growing species. The problem of combined harvesting 
of two ecologically independent species has been studied [10, 13]. The effects of harvesting on the dynamics of 
interacting species have been studied Measterton-Gibbons [14], Chaudhuri et.al. [6-9] with constant harvesting, 
the prey predator model is found to have interesting dynamical behavior including stability, Hopf bifurcation 
and limit cycle [4, 5, 11, 15]. 

The multi species food web models have found to have rich dynamical behavior [16, 18]. S Kumar et. 
al. [17] have investigated the harvesting of predator species predating over two preys. 

In this paper the dynamics of two ecologically independent species which are being harvested have been 
discussed when the dynamics of effort is considered separately. 



II. THE MATHEMATICAL MODEL 

Consider two independent biological species with densities X 1 and X 2 with logistic growth. The 

Mathematical model of two harvesting prey species with effort rate is given by the following system of ordinary 
differential equations: 
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The logistic growth is considered for the two preys. The model does not consider any direct 
competition between the two populations. The constants K iy r iy A , and B, , are model parameters assuming 
only positive values. The effort E is applied to harvest both the species and C is total cost of fishing. The 
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harvesting is proportional to the product of effort E and the fish population density X t .The catch-ability 
coefficients (f L are assumed to be different for the two species. In the model, the third equation considers the 

dynamics of effort E. The constants p x and p 2 are the price of the per unit prey species. The last equation of 
(1) implies that the rate of increase of the effort is proportional to the rate of net economic revenue. The constant 
k is the proportionality constant. 

Let the constant Mo is the reference value of E . Introduce the following dimensionless transformations: 

t = r l T, y t = X i /K i (i = 1,2), x = E/M 0 ,w l =A l q l E 0 /r l , w 2 = B l K l ,w 3 = B 2 K 2 
w 4 = r 2 / r l ,w 5 = A 2 q 2 M /r 10 ,w 6 =kKjM 0 ,w 7 =kK 2 /M 0 ; 

The dimensionless nonlinear system is obtained as: 
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III. EXISTENCE OF EQUILIBRIUM POINTS 

Since 0 < y. < 1; i = 1,2, the underlying non-linear model (2) is bounded and has a unique solution. 
There are at most seven possible equilibrium points of the nonlinear harvesting model: 

E 0 =(0,0,0), £,=(1,0,0), E 2 =(0,1,0), E 3 =(1,1,0), 

£ 4 =(>’,, 0,x*), x* = (1 - y* )(1 + w 2 y* )/ w 1 ; y* 3 =C/(w 1 p 1 w 6 -Cw 2 ) 
E 5 =(0,y* 2 ,x*), x* =(w 4 (1 - y* )(1 + w 3 y 2 ))/w 5 ; y 2 = C/(w 5 p 2 w 7 -Cw 3 ), 

E 6 =[y l ,y 2 , x )• 

Theorem 3.1 The equilibrium point E 4 = ( y, , 0, X ) is feasible only when 
C < w t p t wj(l + w 2 ) 

Theorem 3.2 The equilibrium point E 5 = ^0, y 2 , X j is feasible only when 

C < w 5 p 2 w 7 /(I + w 3 ) 

The proofs of the two theorems are straightforward as 0 < y <1; i =1,2. 

Theorem 3.3 The positive non-zero equilibrium E 6 of nonlinear harvesting model (2) exists provided the 
following conditions are satisfied: 

C < w lPl wJw 2 ; C < w 5 p 2 w 7 / w 3 (5) 

Proof: Proof is given in [6]; the equilibrium points are 

** = (i - y* )(i + w 2 y* + y 2 )/ w i 

(w 5 - w,w 4 )(w 7 p 2 w 5 - Cw 3 ) + Cvv’ w 4 



(3) 



(4) 
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w 5 (w 7 p 2 w 5 - Cw 3 ) + w,w 4 (w 1 /? 1 w 6 - Cw 2 ) 
Cw 5 - (w 5 - WiW 4 )(w 1 p 1 w 6 - Cw 2 ) 



w 5 (xv 7 p 2 w 5 - Cw 3 ) + W l w 4 (w,/7,w 6 - Cw 2 ) 

The positive non-zero biological equilibrium E 6 = ( y x , y 2 , x ) exists provided the conditions (5) are satisfied. 
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It may further be observed that conditions (3) and (4) imply (5), that is if E 4 = ^ y* , 0, Y* ) and 
E 5 — ( 0 , y 2 ,X ^exists then E 6 will also exists. However, the existence of E 6 is possible irrespective of 
E 4 and E 5 provided the condition (5) is satisfied. 



IV. STABILITY ANALYSIS 

The variational matrix about the point E 0 is given by 

10 0 
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0 w, 0 



0 0 -c 

From the above variational matrix, it is seen that there are two unstable manifolds along both X , Y axis and 

one stable manifold along Z axis. Therefore the point E 0 is a saddle point, that is, at very small densities of 

species the effort decreases and tends to zero, while for small efforts the densities of harvesting species will start 
increasing, 

The variational matrices about the axial point E { = (l, 0, 0) and E 2 = (0, 1, 0) are given by 



respectively. 
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From the matrix J x , it is seen that there exists a stable manifold along X axis and an unstable manifolds along 
Z axis. Stable manifold along Y axis exists provided W l p l W 6 — C(1 + W 2 ) < 0 . Observe that this condition 
violates the existence of E 4 = ^ y x , 0, X ) . The point E { is a saddle point. 

Similarly, from the matrix J 2 , it is seen that there exists a stable manifold along Y axis and an unstable 
manifolds along X axis. Stable manifold along Z axis exists provided W x p x W 6 — C(1 + W 2 ) < 0 . This 
condition excludes the existence of equilibrium point E 5 = ^0, y 2 , X* ) .The point E 2 is a saddle point. 

The variational matrix about the point E 3 = is given by 
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Thus, the equilibrium point E 3 = is stable provided the following condition is satisfied: 



W\P\W 6 + w 5 p 2 w 7 
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-C< 0 



( 6 ) 



Theorem 4 . 1 : The equilibrium point E 4 =( Tj*, 0 , X ) is locally asymptotically stable provided 
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( w 2 -l)/2 w 2 < y* < (w, - w 5 )/w 5 < 1 



(7) 



Proof. The variational matrix about the point E 4 =(y 1 , 0,x j is given by 
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The equilibrium point £4 is locally stable if the following conditions are satisfied: 



WjX* > w 2 (l- Tj) 2 ; and y* < (w 1 -w 5 )/w 5 ; 

Substitution for x* and simplification yields the stability conditions as 
((w 2 - 1)/ 2 w 2 ) < >’,* < ( ( w 4 w, - w 5 )/wjW 4 ) < 1 

The equilibrium £4 is unstable when the condition (10) is violated. 

Similarly, the stability conditions for the equilibrium i? 5 are stated in the theorem 4.2. Its proof is omitted. 



Theorem 4.2: The equilibrium point E 5 — (0, y* 2 ,x*^ is locally asymptotically stable provided 

(w 3 -1)/2w 3 <;y* <(w 1 -w 5 )/w 5 <1 (8) 

The equilibrium £5 is unstable when the condition (11) is violated. 

The following theorem gives the conditions for the locally stability of the nonzero positive equilibrium 
Point E 6 = (y*,y*,X*). 



y t ,y 2 ,x ) is locally asymptotically 

stable if the following conditions are satisfied: 

x* > w 2 (l-y*) 2 ; (9) 

w 4 wfx* > w 3 w 5 ( 1 - y*) 2 (10) 

wfx* y ; >(w 1 w 2 y; + w 3 w 5 y 2 )a- y ;) 2 ai) 

w 4 wfx*y* 2 > (w{w 2 y* + w 3 w 5 y 2 )(l - y*) 2 (12) 



Proof: Proof is given in [6]. 

Thus, the positive non-zero biological feasible equilibrium E 6 is locally asymptotically stable if the conditions 
given by (9-12) are satisfied. 

The following theorem gives the conditions for the global stability of the nonzero positive equilibrium point. 

Theorem 4.4 Let the local stability conditions given by (9-12) hold. The positive non-zero biological feasible 
equilibrium E 6 = ( y x , y 2 , X 1 is global stable if the following condition is satisfied: 
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( w 3 - aw 2 w 4 f < 4 aw 4 ((1 + w 2 y* + w 3 y* 2 ) - w 2 )((l + w 2 y\ + w 3 y* 2 ) - w 3 ) 

a = W^ J w 3 p 2 -w£)_ > Q 

w 5 (\v\w (i p x -w 2 C) 

Proof: Proof is given in [6]; 

V. THE OPTIMUM HARVEST POLICY 

The economic rent (net revenue) at any time is given by 
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The present value / of a continuous time stream of revenues is given by the expression 
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Consider the integrand of the present value 
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Here 5 is the instantaneous annual rate of discount. 

Using classical Euler necessary conditions to maximize the positive nonzero equilibrium such that 
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On solving the first equation, we get 
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CVO (24) 

Eliminating 8 , C 1 , C 2 using (26) and (27) we get 

yi~w A y 2 = (\-w A )l2. (25) 

Hence the optimum line (21) and (25) are same in the both cases. Therefore the optimum equilibrium will be the 
intersection of line (25) and first isocline of (2). 

VI. EXISTENCE OF HOPF BIFURCATION 

The characteristic equation of the above variational matrix about E 6 is obtained as [ 6 ]. By the Routh- 

Hurwitz criterion, the positive nonzero equilibrium point is locally asymptotically stable under the conditions 
(9-12). 

The Hopf bifurcation occurs at C = C where the value of a 0 a x - a 2 becomes zero. Substituting a Q a x = a 2 
into (16), we get, 

+ a^/ i 2 + a^y l + a 0 a^ — 0 => (X 2 + a^ )(X + a 0 ) = 0 . 

This gives purely imaginary roots and one real root: 

X± a^ 5 ^,3 — ^ i (26) 

Trans versality condition: - Let the characteristic equation be such that it contains a pair of purely imaginary 
roots X 1 = X l + iX 2 and one real root, say C 1 : 

(X-Xi)((X-X 1 ){X-c x ) = 0. 
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Comparing the coefficients of (31) and (16) we get 
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Differentiating (30) with respect to bifurcation parameter C , setting C= C , and rescaling that 
X , (C ) = 0 .we get [18], 
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Thus the transversality condition is satisfied. So there exists a family of periodic solutions bifurcating from E (> in 
the neighborhood of C* that is, C e(C - S, C +S). 

VII. NUMERICAL SIMULATION 

As the solution of the system is bounded, the long time behavior of the solution is obtained as limit 
cycle, limit point attractor, quasi-periodic. For global dynamic behavior, numerical simulations of the 
underlying non-linear system are carried out. Consider the biological feasible set of parameters as 

w, =2.0, w 2 =1.2, w 3 =1.6, w 4 =1.12 ,w 5 =2.5, w 6 =1.3, w 7 =1.7,^ =0.15, p 2 =0.12, (30) 

The sign change from positive to negative for the expression a^a^ - <X 2 is observed as the values of 

parameter C are varied. The existence of hopf bifurcation is observed in the neighborhood of C = 0.14. 

The variational matrix at the hopf bifurcation point C = 0. 14 is given by 
-0.1725 0.1967 -0.3617“ 



/i = 
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-0.3406 
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The eigenvalues of the matrix J x are obtained as —0.2989, 0.0208 + 0.2443i, 0.0208 - 0.2443i . 
From the eigenvalues the Hopf bifurcation is evident is evident in the neighborhood of C — 0.14. 

The value of the expression Ct 0 a x - a 2 is 0.0030 >0 atC =0.156. Fig. 1 (a) shows the phase 
space trajectories converging to the point (0.42973, 0.36954, 0.59802) for the data set (30) at 
C = 0.156, starting with two different initial conditions. In other words, the nontrivial equilibrium point 
(0.42973, 0.36954, 0.59802) is stable giving the persistence of the system for the given set of 
parameters. Fig 1(b) shows the time series for the data set (30) at C — 0. 156 . 

Figure 2 (a) shows the phase space trajectories and their time series for the data set (30) at C — 0.14. Figure 2 
(b) shows the long time behavior in phase space trajectories and their time series for the same data set (30) 

atC = 0.14. The solution is quasi periodic. The value of the expression a 0 a x ~a 2 is 

-0.0057 < 0, at C = 0.14 
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Fig. 1: For the data set (34); (a) 3D behavior (b) time series 
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Fig. 2 (b). Quasiperiodic behavior 

However, for C = 0.11, there exists a limit cycle for the same set of data. As is evident from the Fig 3. The 
value of the expression CLqCL^ - a 2 is -0.0077 < 0 , at C = 0. 1 1 . 




VIII. CONCLUSIONS 

In this model, separate dynamics of harvesting effort is considered. The local and global persistence of 
the harvested preys has been analyzed. The optimum population level is investigated. The existence of hopf 
bifurcation with respect to total cost of fishing is observed. 
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